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We present a nonperturbative renormalization group solution of the Gell-Mann-Levy cr-model 
which was originally proposed as a phenomenological description of the dynamics of nucleons and 
mesons. In our version of the model the fermions are interpreted as quarks which interact via the a 
and 7r mesons. We derive and numerically solve renormalization group (RG) flow equations to leading 
order in a derivative expansion to study the behavior of the model as it evolves from high to low 
, momentum scales. We develop an expansion in chiral-symmetry-breaking which enables us to track 

■ this symmetry breaking with the evolution of the scale. We use infrared observables to constrain 

the phenomenology allowing predictions of other quantities such as tv — n scattering lengths. The 
. 1 results show improvement over the tree level calculation and are consistent with experiment and 

the results of alternate theoretical approaches such as chiral perturbation theory and lattice gauge 
y-^ \ theory. 
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I. INTRODUCTION 



The complete elucidation of nuclear dynamics in terms of Quantum Chromodynamics (QCD) is still forthcoming. 
This elucidation involves a meeting in the middle of two complementary efforts. On the one hand, using mainly 
, symmetry patterns and the wealth of low energy experimental data, nuclear physicists have been able to construct 
<3\ ' models valid for different regimes and deduce relationships between these models. On the other hand, since the early 
1970s it has become increasingly clear that nuclear processes are the emergent phenomena of the dynamics of quarks 
and gluons (QCD). That these threads can be pursued independently has been demonstrated by the quantitative 
success of nuclear models without any knowledge of QCD and the subsequent widespread acceptance of QCD largely 
on the basis of high energy behavior that appears to have little or no bearing on nuclear phenomena. But in the 
(-H I middle there ought to be some confluence. 

Indeed, though the details are still murky, QCD does dictate that quarks and gluons are confined into hadrons 
within a region of about 1 fcrmi. Thus the QCD degrees of freedom at lower energies are transformed into the 
familiar nuclear degrees of freedom. One can then begin with a model of quarks and mesons at the confinement 



scale and ask what relationship each of the parameters has with QCD parameters. The form of the low energy 
dynamics is strongly constrained by chiral symmetry and its breaking in the hadronic sector. The small departure 
from chiral symmetry is a consequence of the smallness of the up and down quark masses with respect to the QCD 
scale Aqcd ~ 200-300MeV. This can be exploited to construct an effective field theory for low energy QCD or what 
has been termed chiral perturbation theory (\PT). Here the picture is based on a sigma model whose parameters are 
determined phenomenologically. The question arises, however, how to relate the parameters of the effective low energy 
theory to those of short distance QCD. Intermediate to addressing this problem it is a desideratum to compute the 
scale dependence of parameters of simple models which incorporate the chiral symmetry breaking patterns of QCD. 

In the present work we solve the Gell-Mann-Levy cr-model Q] nonperturbatively using Renormalization Group 
(RG) flow equations where projections onto momentum independent couplings at each momentum step are made [0. 
This is equivalent to truncating the derivative expansion to leading order (LO), an approximation sometimes referred 
to as the "local potential approximation". A small quark mass is introduced at the ultraviolet scale to break chiral 
symmetry. The free parameters of the model can be constrained by infrared observables allowing modest predictive 
power. The main result of the calculations, detailed below, are predictions of 7r — n scattering lengths which are 
significant improvements over tree-level calculations and are consistent with experiment as well as with lattice and 
chiral perturbation theory results || . The present paper summarizes many of the results reported in Ref . Q . 

The organization of this paper is as follows. In section [0] we briefly outline the history of the RG approach to field 
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theories. In section [II, our derivation of the RG flow equations for the cr-model is sketched. The numerical technique 
for solving the flow equations is then outlined, our approach to constraining the parameters are discussed, and the 
results presented in section |v|. A summary and some extensions to the model are considered in section VI. Three 
technical appendices include derivations of results used in the text. 



II. REVIEW OF THE RG FLOW EQUATION 



The renormalization group equation upon which this work is based was first derived by Wilson and Kogut || and 
Wegner and Houghton || in two different but equivalent ways. Wilson and Kogut derived their equation using a 
so-called "smooth cutoff" between integrated and unintegrated modes. Their rationale was to avoid nonlocalities 
that naturally arise in position space when one makes sharp divisions in momentum space. The price paid for this 
approach is that unconstrained smoothing functions must be introduced to facilitate the integrations. These smoothing 
functions complicate the equations and make their solution — even after drastic approximations (e.g. the leading order 
(LO) in the derivative expansion.) — impossible analytically. There are also unphysical dependencies on the detailed 
form of the smoothing function. The so-called "sharp-cutoff" method of Wegner and Houghton j(| is simpler in that 
there are no smoothing functions to contend with, but nonlocal effects are now present. Another problem with this 
approach is that there appear to be ambiguities in the form of the RG equations that depend on a variable change at 
a particular point in the derivation p5| . However, at least to leading order (LO) in the derivative expansion, there is 
no ambiguity with the sharp cutoff approach which we employ in this paper. 

One can imagine many different ways to approximate the exact RG equations. A natural approximation scheme 
involves expanding the action in powers of momentum, or — in real space — in powers of derivatives of the fields. After 
integrating over field components corresponding to momenta in the UV momentum shell, the leading order (LO) 
approximation is then obtained by setting all of the remaining components equal to zero except for the uniform 
component. (This is sometimes called the "local potential approximation" (LPA).) The next to leading order (NLO) 
result is obtained by keeping one small momentum and so on. The LO approximation to the full RG equations appears 
first in the work of Nicoll et al. j?J . 

Hasenfratz and Hasenfratz fa], in an early seminal work, showed that the RG equations approximated using a 
method equivalent to LO in the derivative expansion give interesting and nontrivial results. They solve for the RG 
flows of the effective potential in a pure scalar theory and extract critical exponents to compare with calculations 
performed by other means. For d — 3 they find impressive agreement for the critical exponents v and to. Also 
for d = 4 they find no nontrivial fixed point solution which is consistent with the triviality of <f> A field theories. 
Their work showed that even the LO approximation gives a rich quantitative description of the critical properties of 
strongly-coupled scalar field theories. 

There has been a recent resurgence of interest in the RG effective action approach to scalar field theories. Hasenfratz 
and Nager || study the cutoff dependence of the Higgs mass using RG methods. Wetterich and co-workers [^0|~^2|, 
essentially extend the work of Hasenfratz and Hasenfratz to NLO use a smooth cutoff procedure. This allows them 
to include the effects of wavefunction renormalization and compute the critical exponent rj which they find to be in 
rough agreement with other calculations. Morris |l3| ] studies approximations to the exact RG as derived using both 
smooth and sharp cutoffs and in Ref. computes critical exponents to NLO noting that the scheme appears to 
converge. In Ref. (ll| the so-called sharp cutoff ambiguities are treated in detail and the DE of the effective action is 
reviewed in |16|. Afford [lTjl , with an eye toward the electroweak phase transition, computes exponents with a sharp 



cutoff procedure and discusses some of the practical difficulties involved with extending the calculations to NLO. In an 
extensive RG and Monte Carlo analysis of general scalar field theories, Shepard et al. derive "latticized" RG flow 
equations in LO and demonstrate impressive aggreement with Monte Carlo results for a wide variety of cases, broken 
and unbroken phases, three and four dimensions, and 0(1) and 0(2) theories. Some recent related work appears here 
as Refs. 



In addition the inclusion of fermions has been addressed in the literature recently. Maggiorc includes fermions 
via a generalized Yukawa term and derives LO flows for the scalar and generalized Yukawa potentials. He finds no 
evidence that the fixed point structure of the theory is affected by the addition of fermions. In a more comprehensive 
study, Clark and co-workers |2^] derive the exact RG equations for theories of arbitrary field content and derive from 
this equation the LO flows for the scalar and generalized Yukawa potentials which agree with the results of Ref. j|,^2| . 
In two related papers Clark et al. study the issue of computing mass bounds for scalars and fermions in the standard 
model and the stability of fine tuned hierarchies pq| . In related work, Ellwanger and Vergara |2(| use RG flow 
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equations for generalized NJL models to study the Higgs top quark system to leading order in 1/N C . Other work 
pertaining to the Higgs top system appears as Ref. |27j| . 

In an effort paralleling this work Jungnickel and Wetterich have applied renormalization group methods to 
the generalized cr-model. In one sense their work is more ambitious than the present work. They have considered 
the general SU(N) sector with all allowable functions of the field invariants and included via a smooth cut-off the 
NLO quantities which allow treatment of the anomalous dimension. But the price paid for such generality is an 
abundance of free parameters, some non-physical, which are difficult to constrain. In contrast, the present work 
focuses exclusively on the SU(2) sector which requires only one field invariant, uses a sharp cut-off approach as in Ref. 
H truncated at leading order in the derivative expansion. While this precludes any treatment of the strange-quark 
sector or the anomalous dimension, there are fewer parameters which can be well constrained by experimental data 
allowing modest predictive power for the theory. 

There is now ample evidence in the literature that the leading order approxiamtion to the RG flow of a field theory 
incorporates correctly many of the intricacies of strong-coupled quantum field theories. The aim of the present work 
is to apply this powerful method to the the Gell-Mann-Levy cr-model in an attempt to understand better low energy 
nuclear phenomenology. 



III. DERIVATION OF THE FLOW EQUATIONS 



A. General 



We first present a derivation of an exact RG equation for an arbitrarily complicated quantum field theory described 
by the action, S( A )[$], regulated at some large momentum cutoff A. The field content of the theory can in principle 
be anything at all, $ = {0, ip, A^, F^, . . .} i.e. fields described by complex scalars, spinors, vectors and/or tensors, 
though we will only be concerned with theories containing scalars and spinors. We define the effective action at the 
momentum scale A — AA, ( 5 , ( A ^ AA ) a g 5 



(1) 



shell 



where the subscript u shell!' indicates that only the Fourier components of the fields {& q } with momenta in the shell, 
A — AA < |g| < A are integrated. Thus the notation for the measure means, 



x>$x>$ = Y[ d<f> qi Yl 



<I2 ■ 



(2) 
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for A — AA < |c/i|,|c/2| < A. $ is the "conjugate" of <& {e.g. complex conjugate for the complex scalar field 4>, 
Dirac adjoint for the Dirac spinor ip, Hermitian conjugate for matrix fields etc.). The actions 5 , ( A ~ AA ) and 5^ A ' are 
considered equivalent in the sense that for \q\ <C A they each give the same n-point functions. We now decompose 
&(x) into uniform and nonuniform pieces, 



$0) = Q a +Cl(x) 



(3) 



Expanding S^' in a functional Taylor series about the uniform field components gives, 



^ A )[$,$] = 5^[n ,n ] 



q^O 



5S 



sn 



^ ss 
+ "<j — 

0\l a 



+ ' n <H 

- jn + Tij + Tlmvl. 



5 2 S 



5Q qi Sfl q2 



(4) 
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The subscript zero means that all modes, Q q , with q in the shell are set to zero; and the primed summation symbol 
indicates that only momenta in the shell are included in the sum. For the purposes of deriving differential equations 
with respect to the independent variable A, we will eventually take the limit AA — > 0, thus we can truncate the series 
after the quadratic term without approximation since all higher order contributions to the RG flows will be at least 
0(AA 2 ). In the second equality several definitions have been made: 



/ ss 
I sn„ 



J = 



sn n lo 

SS I 
<50 00 lo 



\ 



V 



(5) 
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and its conjugate are the "generalized source" column vectors and the supermatrix of second derivatives is defined as 
(see Appendix A for a review of the supermatrix formalism) , 



M 



£ A 
A T 



(G) 



Also, where momentum subscripts are not present, matrix multiplication over the momentum indices is implied, e.g., 



SS I SS I 

sn qi lo sn q2 lo 



V 



V — 



(7) 



Now substituting Eq.(^) into Eq.(|l|) gives, 

_(j(A-AA) 



1° / Vlivn e -mm + jn+nj) 



S W \o {e JM-iJ sdet -l Mh 



Irrelevant constant factors have been dropped wherever they appear. Using the identity, 



—str In M 



sdet _1 M 

the exact RG equation for this generalized system is easily obtained as, 



S (A-AA) = S (A) 



str In M - JM- X J + C(AA 2 



(8) 
(9) 

(10) 



This equation relates the action at momentum scale A — AA to the action at momentum scale A. We will not convert 



23[ | do) until we consider its form for a particular 
The "str In" term is generally referred to as the 



it directly into a functional differential equation (as Clark et 
action. A terminology has grown up around equations such as ([10 
"loop" term and the "3M~ X J" term is referred to as the "tree" term. This is because the contribution from the 
JM~ X J term is present in mean field theory whereas the contribution from the str In term includes effects from loop 
integrations p9| ]. 

For the purposes of treating the tr-model we consider only real scalar and spinor fields (i.e. $ = {4>, i[>}). If the flow 
equations are truncated to LO in the derivative expansion, many simplifications ensue. The source column vector 
becomes, 



J= (i 



V i J 



(11) 
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where, 



SS 



vifjq 
SS 



(12) 



The matrix of second derivatives is, 



M 



£ A 
A T 



S 2 S 



s 2 s 



5tj>_ Stji q 



(13) 



Now we can rewrite the superdeterminant appearing in Eq. (|l^) using some tricks from Appendix |X| (compare with 

Eq.rtAn])), 



sdct _1 M = (dct" 1/2 7V)(det^") 



(14) 



(where N = S — AT 1 A) which gives, 



strlnM = -trlniV - trlnj". 
2 



(15) 



(The 1/2 coming from the counting in the momentum sums due to the real <j).) Now Eq.fllCf) becomes 



5 (A-AA) = g(A) 



-trlniV - trhxF - JM~ X J. 



(16) 



Restricted to LO in derivatives the last ("tree") term in this equation may be neglected. One way to see this 
diagrammatically is to note that any particular contribution to the tree term in Eq. (|10j) or (|16|) will be comprised of 
external legs connected by a single internal propagator with momentum in the shell (see Fig. 111). Since to LO we set 
the momentum of all external legs to zero, these diagrams cannot conserve momentum at their vertices and therefore 
vanish. To make the argument analytically we write a particular term of S' q , 



SS 



-V'" ^2 0?i^? 2 ^3i+g2+«',O H 



2! 



(17) 
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where q' is in the shell. When, in going to LO we set all field components, (f> q , with nonzero momentum to zero, it's 
clear that this term and all higher order terms will vanish. Thus we write Eq.(|l6|) as, 



S (A-AA) = S (A) 



--trlniV + trlii.? 7 . 
» 2 



(18) 



This is the form of the RG flow equation that we use for our calculation in the present work. 



B. CT-model Flow Equations 

We now apply the results of the previous section to derive the flow equations for the cr-model. The incorporation 
of chiral symmetry breaking brings a substantial price in the complexity of the algebra. Many of the tedious details 
are relegated to the Appendices. Ref. 0] includes many details not related in this paper. 



The Euclidean action for the cr-model is: 
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= J d d x 

A few definitions to compactify the notation have been made 



V (A} (p, a) + \{d^ a f + V# + U (A) (f, a, 0)]if, 



(19) 





"(5)- 






pa 


-( 1 . 






g 


= r>° = 


p 2 


= ee f = ( 



a = 0,i; i = l,2,...,iV 



(7 + 17 T ■ 7T, 



(20) 



where the last term is the only field invariant when specifying r a only to the N — 4 or 0(4) case. We will often use 
the more general notation but will only be considering 0(4) in this work. Recall that SU(2) x SU(2) rotations of 
the fermions induce 0(4) rotations in the boson sector. The scalar sector of the c-model then is a four component 
field. The potentials V^(p, a) and lf( A '(p, a, Q) include terms proportional to powers of the zeroth component of the 
scalar field, a, which will break chiral symmetry. To keep track of the symmetry breaking as the theory flows from 
high to low momentum we expand the potentials in the form, 



V(p,a) = V (p) + vVi{p) + —V 2 (p) + 
U (p, cr, g) = m(p, a) + gg(p, a) 



(21) 



where m{p,a) and g{p,a) are expanded similarly to V(p,a). Each of the functions Vk, mu and gu for k = 0,1,2 is 
SU{2) x SU{2) symmetric since they only depend on p. Chiral symmetry is broken by Vi, V%, gi, 32 and m(p,a). 

The LO flow of the action Eq.([l9|) is given by Eq.(|lJ) rewritten as, 



S (A-AA) = S (A) 



-trlnN + trhxF. 



(22) 



where N has components N' 



ah 



N 



S A 
A T 



ab 



A^-'A" 



(23) 



s 2 s 



"-q f q 

S 2 S 



s 2 s 





S 2 S 



5ip_ q 6ip q 



0/ 



The extra indices now refer to the components of the scalar field. We Taylor expand the potentials, 
V w (^ a ) = V {A) (#}) + V' wa (^)cp a + V'( A )° 6 ($J)^V + ■■■ 



(and similarly for U). With normalization 



d d x e lq - x = (Vo\)5 q , 



(24) 



(25) 



we get, 



Vol 



+ E \l<t> a - q [5 a \ 2 + v"^ a \^) + %u"^ ab (r )M <t> h q 

-H£ q [it + uW(4%)]ii% 

+u'^ a (r )(^ a - q ^+^- q ^ ) 



(26) 



G 



From this expression we can compute all the matrices in Eq.(|23), 



Sab cab 2 , r /' ab / ia\ 

9Tl' a Tl' b 

n ab = u ab 



H + u 



N 

Vof 



(27) 



With these definitions Eq.(|22j) becomes, 

C(A-AA) 

Vol 



y( A-AA) ((/)g)+ ^ oC/ (A-AA) (0 a^ o 

^ (A) (^)+?o^ (A) (^)^o 



2Vol 



-tr In 



S a6 + V « a >o 



— tr ln^+[/( A >(0 



(28) 



modifications. First, we can write, 



We can write Eq.(|28[) as two coupled equations for the flow of V and U. Before doing so, we make a number of 

H + U _ Eft 

(29) 



1 

ii + U ~ q 2 + UW ^ IV 



D F = q 2 + UU^ = q 2 +m 2 + 2<rmg + p 2 g 2 , 



where we've dropped the term proportional to ^ since it will vanish at LO when wc take angle averages. Second, we 
can write the trace over the inverse fermion propagator term as, 



tr ln(^ + U) = -nfCd h\D F 



(30) 



with Cd = 2 d / 2 {2^ d x )/ 2 ) for d even (odd), using the identity + U\ 2 = q 2 + U 2 ; the njCd ^2 q factor comes from the 
trace over the flavor and Dirac indices. Also we can write, 



tr In 



= In dct £ ab + tr [(£° c ) _ Vo n M 



(31) 



where the trace in the last term is only over the ab indices. Putting all this together, we can write the V and U 
equations from Eq.(ES) as, 



^A-AA) = y(A) _ _i_ ^(i n det £ - n f c d In D F ) 



jj(a-aa) = ^(A) _J_y tI s -i . n 

2Vol ^ 



(32) 



Now taking the limit AA — > we have, 



dA 



-^A rf (In dct £ - n f c d In D F ) 

-^A d tr (S- 1 • H), 
2 v ; 



(33) 
(34) 



7 



where the determinant and trace are only over the ab indices. In Eqs.(|33|) and (|3~i|), 

- v" ab 
2U' a Wu' b 



^06 ^ab ^2 I -y ab 



n ab = u" ab - 

Dp 

Dp = A 2 + m 2 + p 2 + 2amg. (35) 

The matrices V ab , U ab , and U a are worked out in terms of derivatives with respect to p in Appendix B. One 
glance at these expressions is adequate to impress the reader of the proliferation of algebraic complication in the 
extension of the simple Yukawa coupled fermions jl 22 25| to the case of broken chiral symmetry. The matrices can 



be simplified by a similarity transformation, however, since the determinant and the trace are invariant with respect 
to such transformations. As discussed in Appendix B, when similarity transformed, £ has only six nonzero elements, 



£ = S ■ £ • S T = I "i" " " I, (36) 



Eoo 


£01 








£10 


En 














£22 














E33 



but ft still has every element nonzero although each of the elements is somewhat simpler. Notice, however, that since 
only the combination tr£ _1 • ft appears in the flow equations, only six elements of ft will contribute so that, 

tr S'- 1 • ft = (S 00 )- 1 ^ 00 + (S 10 )- 1 ^ 01 + (E 01 )- 1 ^ 00 

+ (S 11 )- 1 r> 11 + (E 22 )- 1 ^ 22 + (E 33 )- 1 ^ 33 , (37) 

where, e.g., (E 00 )^ 1 means the 00 element of £ _1 . 

We can now consider the "chiral limit" where we set all chiral breaking terms (m(p,a), Vi, V2, gi and g%) to zero. 
Then, 

E afc — > (A 2 + Vo(p)/p)S ab + (V "(p) - Vo/ P )^, (38) 

P 

which becomes, after a similarity transformation, 

E a "' =fn 7 n n I > ( 39 ) 



where A = A 2 + Vq" and A = A 2 + V Q /p. Thus, 



and so in the chiral limit only the diagonal elements of ft will contribute to the trace term. These can be shown to 
be, 

o00 r>\ " 1 2(70 / 2 1 n 2 2 1 2 '2m 

p Of 

n 11 = n 22 = n 33 = g(5a + M). r 4 i) 

Thus from Eqs. (|33|) and (|34|) the 0(4) flow equations in the chiral limit are, 



A 














A 




























A 



AI 1 














AT 1 














A7 1 














At 
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A 



A 



dA 
dgp 
OA 



^A rf (ln L> + 31n - n f c d In D F ) 



-^a4± 

3 



.?;+^-|f(.9o 2 +2 5 oV+A; 2 ) 

P U F 



(42) 



(43) 



For O(N) symmetric theories there would be an A^ — 1 in place of the 3s in front of the InD^ and the (l/D w )(- ■ ■) 
terms. In Ref. JUpHJlS}] the flow equation for an O(N) scalar-only field theory is derived and is equivalent to Eq.(^) 
with rif — and A" = 4. Also Eqs.(|4|H) represent a generalization of the flow equations for the simple Yukawa 
coupled fermion models derived in Refs. ]4|, |2"2| , |23| ] to a model with both scalar and pseudoscalar bosons coupled to 
fermions. 



IV. NUMERICAL RESULTS 



We now discuss the numerical solution of the LO cr-model RG flow equations, ( p3|) and (pi|). As mentioned in 
Appendix B, these equations is quite lengthy. Since we keep only 0(a) (first order in the chiral symmetry breaking 
parameter) in Eq. (plf), there are six coupled flow equations, which for 



det £ = F (p) + vF 1 {p) 
tr E'- 1 • ft = Am Q (p) + QAg Q {p) + a 



A mi (p) + gA gi (p) 



(44) 



take the form: 



OA 


-^±A d (\nF {p)-n f c d In D F 


dV 1 
~dA = 


Ad^F^p) 
2 F Q (p) 


dm 
dA 




drm = 

dA 


-^A d A mi (p) 


dA 


-^A d Ag (p) 


A^i = 
dA 


^A d A 9l (p). 



(45) 



Symbolic computing (we used Mathematica) greatly facilitates the determination of the functions Fq, Fi, Amo etc. 
which are too lengthy to be reproduced here [JjJ]. (Agi(p) and Am\(p), for instance, require several pages of output!). 
Once these functions are determined, however, the numerical solution of Eqs.fljq) proceeds exactly as in the scalar 
only case as reported, e.g., in Refs f§||l8|| except with six functions instead of two. There are a number of nonphysical 
subtleties associated with the numerical solution of Eq.(45) (e.g. the domain limits, the algorithm used to solve the 
equations, and the number of terms retained in the polynomial fits , etc.) These nonphysical sensitivities yield a small 
spread of output values for a given input. At present we have three independently constructed codes to solve the 
equations which give essentially the same results which gives us some confidence in the solutions presented here. 

We set our "initial condition" by specifying the value of the parameters at the UV scale: 



)(p,a) = \pV + \\ p 4 
U^(p,a,g)=m° Q +g g. 



(46) 



Then we expand in powers of p, 
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M 

1i 



^ (A) (p)=E^f ) (A) 



2T" 

i=l 

M 

,( A ), 

//I < ' i I . \ I — 7 

2i 



-r j (p)=d o) (A)+5:^r(A)P 2i 



5i A) (p) = yi 0) (A) + E^f ) (A)P 2l I (47) 
for fc = 0, 1, 2 (see Eq.(^T|)). Thus at the UV scale we have 



^ 2) (A ) 


= Mo 


4 4) (Ao) 


= Aq 


^ 0) (A ) 


= m° 


^ 0) (A ) 


= .9o- 



(48) 
(49) 

with all higher order coefficients for k = 1, 2 set to zero. Just as in the previous section we perform a fit of the functions 
Vfe, rrik, and to a power series in p at each A step. Also, since we're interested in spontaneous symmetry breaking, 
we will have the parameter f n =< <\> > vac ~ &, which sets the scale of the symmetry breaking. Thus we consider the 
set {m^, /„•, /ip, Ao, m q , go, Aq, A/#} the input parameters to the model. This approach is philosophically different 
from that taken by Jungnickel and Wetterich [O] where /„• was the quantity they wished to predict. 

The basic approach adopted in the present work for constraining the phenomenology is as follows. Since there are 
as yet no experimental or reliable theoretical constraints at the UV scale, we are forced to look to observables at 
the IR scale to determine the free parameters of the model. Perhaps the two parameters in the model most tightly 
constrained by low energy data are m n and f n . Thus we will tune other parameters of the model to get m n and f n at 
their experimental values at the infrared scale. This still leaves four out of the six parameters unconstrained, and we 
must decide which out of {/1q, Ao, mj, go, Ao, Am} to fix and which to tune to m v and f^. Results for scalar-only 
calculations |8"Hl§ll as well as for Yukawa coupled fermions j|,^2],^3) indicate that the two parameters Ao and /ig are 
not truly independent, i.e. we can tune to particular values of IR parameters with large number of values for Ao and 
Po so long as Ao > and large. Thus we fix Ao = 10 and use po to tune the IR parameters. From our experience 
choosing a different value of Ao will only require a re-tuning of po to obtain equivalent results. Next there is the bare 
current quark mass, m q . Experimentally, there is about ±5MeV spread of the values for m u and [^TJ. We will 

tune mj to get the average of the means of these values at the IR scale, i.e. m\^ IR ^ = m q = \ {m u + rn c i) — 7.5MeV. 
We also have the two cut-off scales Ao and Am, only one of which need be tuned with the other fixed since the RG 
equations only care about the ratio. We fix A/^ = m T[ ± ~ 140MeV and will allow Ao to be adjusted. This leaves only 
go- We know that the quark- meson coupling is strong at these energy scales and the Goldberger-Treiman relation 
for the constituent quark model gives g ~ M nuc i eon /3f w ~ 3.366. But this is not an approrpiate constraint in our 
zero-density model since we have no nucleons. We therefore perform three fits, one each for g >,:_,< 3.366 and 
compare the results. To summarize, for each of these values of go we tune {po-, m q , Ao} to {m^, /„■, m q }. Once all 
the parameters are fixed, the model predicts other IR quantities such as m a (sigma mass), \^ v t, \3-pt (the sigma 
four- and three-point couplings), g, a[j, and Oq. Not all these quantities have experimental determinations however; 
the quantities best determined by experiment are the irir scattering lengths aj] and a§ discussed in more detail later 
(see Appendix C for background). 

The results of the three fits for g >, ~, < 3.366 are displayed in Table |. The middle column is the result for fitting 
go to give g ~ 3.36; the right and left columns are the results for arbitrarily choosing go = 2.500 and 3.100 respectively. 
The last three rows (in the top section) contain the actual fits to m n , f^, and m q . The point was to get (m n , m q ) 
~ (140MeV, 92.4MeV, 7.5MeV) & The f and || indicate that m° and (fx , A ) were used to fix m q and {m„,f„). 



The values of m q and m q all fall within the uncertainty quoted in the particle data book 1 30 1 . Rows 9 through 15 
represent some of the predictions of the calculation. Since the sigma is less a "particle" and more a broad resonance, 
the values for m a , \^ pt , and \s- p t are hard to compare quantitatively with experiment. The values for the scalar 
density, < tpip > are close to other calculated values of < ipip >~ —[(240 ± 25)MeV] 3 pi]] . In addition the adjusted 
values of Ao ~ 940MeV are in the range expected. 

Fig. |^ shows the boson potentials as a function of p. As expected for weakly broken chiral symmetry, the first order 
term, V\ is just a small correction. One can see clearly that the minimum is at p =< a > vac = ~ 93 McV [B0[ . 
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Figs. |^ and Q display the other functions computed in the model. At p ~ 93 MeV the values of m and g are just the 
IR values quoted in Table [i]. The relatively small contribution from g\ and m justify ex post facto our expansion in 
the (small) chiral symmetry breaking parameter a (Eq. (pi]) ) . 



V. tttt SCATTERING LENGTHS 



Another set of predictions from the RG solution of the er-model come in the form of the parametrization of low 
energy tttt scattering. The expansion of the real part of the partial-wave amplitude can be written as 



Re,4f( S ) = 327r (JL) (al + bl^ + ■■■), (50) 



2 



where / = 0,1,2 denotes the isospin channel and I is the partial wave index. (See, e.g., Refs. [B2j section VI-4 and 
p3[ . Appendix C contains a review of the perturbative calculation and the connection to our model.) So for low 
energy (q 2 -C m 2 ) scattering af and 6f will be the most relevant quantities to study. Table || displays a comparison 
of our three fits and a number of different calculations and experimental values for and Oq in dimensions of inverse 
pion mass. Also, the quantity 2a[J — 5ag is included since for s-wave scattering it provides a constraint 0]. Each of the 
three fits give results that are consistent with experiment and with x?T and lattice QCD calculations. We discuss 
how the 7T7r scattering lengths are computed in our model in Appendix C. 



VI. SUMMARY AND POSSIBLE FUTURE EXTENSIONS OF THE MODEL 



We present a nonperturbative solution of the cr-model using the a sharp-cutoff RG equation truncated to LO in 
the derivative expansion. The model mimics the chiral symmetry of QCD. An important feature of our approach is 
that we can straight-forwardly track the chiral-symmetry-breaking. Including a small quark mass at the UV scale, 
we can follow the chiral symmetry breaking in the numerical solution as the scale is lowered. The parameters of the 
model can be constrained by low-energy data allowing predictions of, for example, tt — n scattering lengths. The 
values for tttt scattering lengths obtained in this way show an improvement over the perturbative calculation || and 
are essentially consistent with experiment and other nonperturbative calculations |35|]36|| . It might be added that 
these calculations were performed at a substantially lower computational cost than those of Refs. j3^j3(|. This work 
contrasts with that of Ref. [0 in that we confine our attention to the SU(2)xSU(2) sector, use a sharp cut-off approach 
truncated at LO, and allow for explict treatment of the chiral symmetry breaking via an expansion in the vacuum 
expectation value of the cr-field. We also adopt a different philosophy as regards the phenomenology constraining the 
free parameters. 

There are a number of possible extensions to the present calculation. Perhaps the easiest is the extension to 0(cr 2 ). 
Indeed much of the analytical work has already been done with the expected result that, for small current quark 
masses, the second order potentials will be small corrections to the first order, testifying to the convergence of the 
method. We have performed 0{a 2 ) calculations with only the bosonic potentials flowing {i.e. fixing m(p) = m® and 
g(p) — go for all A) and confirmed that V2 is small compared to V\. The results for the other calculated parameters 
in the model are not expected to change substantially at 0(a 2 ) since m° q is small. 

Another possible extension would be to allow for finite density which would permit an analysis of nuclear phenomena. 
Crudely, at finite density, there would be a momentum scale, kp, below which the fermion parts of the flow equations 
would cease to contribute due to Pauli blocking, while the boson parts would still contribute to the flow. Such crude 
calculations using this scheme show the qualitative restoration of chiral symmetry as kp increases. A proper handling 
of finite density with the RG, however, for relativistic field theories requires much more care so we defer offering any 
conclusions at the present time. 

Yet another extension is to incorporate strangeness into the model by extending to SU(3) x SU (3). This requires the 
introduction of two more field invariants as treated in Ref. fl2|| . In this case, since m s ~ Aqcd, chiral SU(3) x SU(3) 
is strongly broken and the expansion in the chiral symmetry breaking field may not converge rapidly enough and some 
other approach to following the chiral symmetry breaking may be required. 
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APPENDIX A: SUPERMATRIX FORMALISM 

In this appendix we review the supermatrix formalism. A "supermatrix" is a square matrix of the form, 

tut - ( m bb M BF \ f An 

m -\m fb M FF )> (A1) 

where the square sub-matrices Mbb and Mpp are both even elements of the Grassmann algebra while Mbf and Mfb 
are odd elements of the Grassmann algebra. 

(Grassmann or anticommuting variables allow the incorporation of Fermi statistics into the path integral formalism 
of quantum field theory. An arbitrary function of a Grassmann variable can be Taylor expanded as 

f{9) = a + [30 (A2) 

where a and /3 are normal numbers obeys the Grassmann algebra, 

{0, 9} = 69 + 90 = => 9 2 = (A3) 

which is why the Taylor series terminates. For more details see , e.g., Ref. [[37| p. 214-219.) 
The supertrace is defined as 

strM = trM B B - tiM F F ; ( A4) 

so that the familiar commutative property of the normal trace still holds for the supertrace. The superdeterminant is 
defined as 

sdetM = e strlnM , (A5) 

which preserves the familiar property of determinants, det M N = detMdetiV. 
Now consider the decomposition of M: 

— (is o(i m ^ f ) m 

where Nff = Mpp — MfbM^ b Mbf- Now it's easy to show 

sdet^^ f)=tetM B B (A7) 

and, 

sdet ( I Mb }^ BF ) = detNpp (A8) 



F 



so that, 

sdetM = detMsijdet^iVpp,. (A9) 
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Similarly we could have chosen the decomposition 



M = ( o 1 FF (m fb M PF > ( A1 °) 



_ (N BB M BF M FF \ ( 1 
\ I ) \M F 

where N BB — M BB — M BF Mp F M FB . Then we would be led to 



sdetM = det7V BB det 1 M FF . (All) 



APPENDIX B: DERIVATIVES OF THE CHIRAL FUNCTIONS 

In the derivation of the flow equations for the cr-modcl, derivatives with respect to <j> a of functions of p, a, and 
Q = T a (j) a are taken. In this appendix we work out what these derivatives are. Consider the potential 

U(p,a,g)=m(p,a)+gg(p,a) (Bl) 



where, 



and similarly for g. Now 



2 

m(p, cr) = TOo(p) + <rm 1 (p) + —m 2 (p) H (B2) 



^L = U a = m a + Gg' a + T a g 

8 2 U „ // / 

U ab = m ab + Qg ab + T a g b + g a T b . (B3) 



d<j> a d<j> b 

We need to work out the derivatives of functions of p and a. Consider f(p, a), f = m or g: 

9£_ = /a = f{0 ya + afW a + ^ f{ 2 Y a + 5 aO (/ (l) + ^(2)) 



= ^/ +(5 a0 (/ (l) +(7/ (2) ) 
P 

f"ab = f(0)"ab + a f(l)"ab + °f_ f(2)"ab + £aO(y(l)'& + a fW'bj 

+(/ (1) ' a + <7/ (2) ' a )£ 06 + 5 a0 5 0fo / < 2 ) 



p p \ p 



+ S^£±£^^ f{1 y +af{2 y^ +sa0s0bf(2) (B4) 



In the last of each of the above equalities we've used 



h a = — h (p) 
P 



h 6 a 6" h 



P P 2 V P 



h I*' 1 - ■ —^[h - ). (Br,) 
where h = h(p), h — mM or <?W for i = 0, 1, 2. Also 
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/(p) = / (0) '+<r/ (1) ' + y/ (2) ' 

f"(p) = f (0) " + *f w " + Yf (2) "- 



(B6) 



With these derivatives, we can compute the matrices in section III B 
S Qb = A 2 S ab + v" ab { Pl a) 

A 2 + — ) S ab + ! I ' 



V \cf) 



,4* + s^ + s, ^ + ^ + y(W 



p 2 p 



(B7) 



and 



Sl Qb = U" ab (p, a, G) - j^U^p, a, g)rf (p, a, g)U b { Pl a, Q) 
Dp 



u 



§ ab 



u 



u ~ 



if 2£/ t ,,\ 6 a 6 b 



D F 



-U 



2W > ~\S aQ <j) b + (/) a 5 0b 
'D F ~ ) p 



2Z/t / c/) a S ab 2C/t , 5 a °4> b 

—f—Qg m h -^—ym g 

D F p Dp p 

Dp J \ Dp J p 

~ q _ ?Elu] ( r a s ob + s a0 r b ) - ^r a u^r b 



i°r" 



ab 



Dp 



Dp 



s <5 a( V b + a <^ oh 
' 1 1 



'50" 



+ "<A5 1" "<5 



+fW a0 <5 0b + n. 



r a <A b 



b a r b 



50- 



a X oh , s a0 T b )-T a Q,aaT h , 



P 



(B8) 



where 



X = X {1) +aX {2 \ X = V,U,m,g 

I I I 

U = m +Gg ■ 



(B9) 



Note that the Ss are all functions of p and a and the Sis are all functions of p, a and Q. The ordering in all the terms 
containing Qs is nontrivial since the Pauli matrices don't commute with each other. 

Now using the similarity transformation matrix (defining it 2 = ix\ + it 2 and 7r 2 2 = it 2 + 7r|), 



t I, 



7Tl 



7T 2 



7T3 \ 

<T"7T-3 



V o o - P | p|/ 



(BIO) 



(obtained by multiplying rotation matrices in 4-space — see Appendix D in Rcf. |4j]) we can compute, 

S' = S • S • S T 

n' = s • n • s T , 



(Bll) 
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this amounts to performing similarity transformations on each of the tensors 

b r£ rs^ 00b 

o , 9 , , ,00, 

p p p 



T a rh h A- d> a V b 

9 +9 — , (r a s ob + s a0 r b ), ft 6 , (Bi2) 



as can be read off from Eqs.(B7) and (JB8|). After performing the similarity transformation X has the form 



X = 



/ ^00 


Eoi 








Eio 


Eii 














E22 





V o 








E33 



(B13) 



since E = E we can write 

det x' = (E 00 E n - (E 01 ) 2 )S 22 S 33 
tr S'- 1 ■ ft = (E 00 )^ 00 + (E 10 )- 1 ^ 01 + (E 01 )- 1 ^ 00 

+ (E 11 )- 1 ft 11 + (E 22 )- 1 ^ 22 + (E 33 )^ 33 . (B14) 



Thus only 12 elements — 6 from X and 6 from f2 — need be computed to determine the flow equations for the a- 
model. As is probably clear this substantially reduces the complexity of the expressions but they they are still quite 
complicated. The derivation is facilitated by using Mathematica to compute the E and CI functions in Eqs.(B7) and 
(B8). Then Eqs.(B14) can be written in terms of these expressions. 



APPENDIX C: tttv SCATTERING LENGTHS IN THE RG MODEL 



In this appendix we will sketch the calculation of tttt scattering lengths. Before addressing the calculation in our 
model we first discuss the calculation using the perturbative a- model. This was first done by Weinberg B; useful 
reviews appear here as Refs. [||J33]. 



Differential cross sections in field theory are computed by squaring the "invariant amplitude" , A which is usually 
computed to a given order in perturbation theory using diagrammatic techniques (for details on our conventions see 
e.g. Ref. || Appendix A-3): 



Consider the amplitude for processes involving two pions in and two pions out. 

As indicated, it will depend on the isospin channel / = 0, 1, 2, the Mandelstam momentum variables s — (pi +P2) 2 , 
t = (pi +P3) 2 , u = (pi + pa) 2 and the isospin indices of each of the pions, i,j, fc, I. (In the center-of-mass frame these 
become s = 4(m 2 + k 2 ), t = —2(1 — cos8)k 2 , and u = —2(1 + cos8)k 2 , where k is the 3-momentum of the incident 
pion and 9 is the angle between the incident pion and the out going pion.) We can expand the amplitude in partial 
waves 

00 

A z (s,t,u) =^(2Z + l)Af(s,t,if)P;(cos0), (C2) 
z=o 



where Pi(cosd) are the Legendre polynomials. Amplitudes for each of the isospin channels are not independent, 
however; since all the particles are bosons Aj. kl (s,t,u) is totally symmetric in all indices. This can be exploited to 
show that there is really only one independent amplitude Ai(s,t,u). Each of the amplitudes A\ can be written in 
terms of Ai(s, i, u): 
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Al(s,t,u) = 3Ai(s,t,u) + Ai(t,s,u) + A t (u,t,s) (C3) 
A\{s,t,u) = Ai(t,s,u) - A t (u,t,s) 
Af(s, t, u) = A t (t, s, u) + Ai(u, t, s). 

So we need only compute Ai(s, t, u). 

At tree level, or lowest order in the coupling constant A, the calculation is quite simple. The Feynman rules can be 
read off of the boson potential, 

V(p) = \p 2 p 2 + \\p A = ■■■- Xcj mm an 2 - ^tt 4 + ■ • ■ . (C4) 
We consider the process tt + it~ — > tt q tt° with 7r = tt 3 , tt± = "^( 7I ' 1 ± tt 2 ), then we have 

V(p) = Xcr min a(7r 2 + 27r + 7r_) - ^(tTq + 4ti-2 71-2 + 47t a 7r+7r _) + . . . . (C5) 



The amplitude can be written as 



Ai(s,t,u) = -2i\ + 4(-i\a min ) % — (C6) 

s — m% 



f 2 J\ m%-s ) f 2 
where we used ml <C m 2 and s = Ami (at threshold) in the last approximation. Also we used the relations 

Cmin = — fir (C7) 

A - (C8) 
P? = \{ml-Zml) (C9) 

to replace (a m i n , A, p 2 ) with the observable (m a , m^, f v ). From Eq. flC4|) we can compute the amplitudes for the 
isospin channels, 

A$(s, t, u) = 3A(s, t, u) + A(t, s, u) + A(u, t, s) (CIO) 
A\(s, t, u) — A(t, s, u) — A(u, t, s) 
Ao(s, t, u) = A(t, s, u) + A(u, t, s). 

using Eq.(|C7|), s + 1 + u = Am 2 , s = A(m 2 + k 2 ) 

4 = ^ = ^ + 8^ (Gil) 

, -i i — u 4 rn?. P 
A 1 - - 71 



1 " fl 37T /* TOj 

2 t + u - 2ml _ 4 to 2 ^ 

/* /2 7T/2TO2- 



The "scattering length" a\ and "slope parameter" are defined by 



Re A\ (s) = 32tt 



mi 



Of- 



(C12) 
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which gives the tree level values [BE 



J_mi_ r 
32tt fix 2 " 
1 ml 



1 m 
167T fir 2 



24tt /tt 2 

2 



0.16 



0.0 



-0.045 



&2 



32tt /tt 2 



0.18 



(C13) 



6} = 



bl 



1 



167T /7T 2 



-0.09 



The calculation in our model proceeds similarly. Since we compute the potential V(p) we must relate the tree 
diagrams in the invariant amplitude to this potential. In the leading order (LO) approximation, all the momenta on 
the external legs of the diagrams are zero and therefore we can only compute the s-wave (/ = 0) scattering lengths. 

We equate the vertex between two pions and one sigma with the third derivative of the potential and the 4-pion 
vertex with the fourth derivative. 



Computation of these is straightforward using Eq. ( B4 ) as a starting point 



d 3 V 



d 4 V 



7? = 



S* 3 



V V 

frr f% 



(C14) 



dir 1 diri dir k dir 1 



s=0 = {S^S kl + 5 %k 5 jl + S il S 3 ' k ) ( ^- - ^ 

J 7T J TT 



where 



V = V (p) + aV 1 (p) + -V 2 (p) 



(C15) 



and similarly for V . Defining F(p) = X and G(p) — F(p)/f„ + F 2 (p) /(p — m 2 ) we can construct the invariant 

amplitude, 



A(s, t, u) = 6 l] 6 kt G{s) + 5 lk 5 ]l G{t) + 5 ll V k G{u) 



(C16) 



At threshold s = Am 2 , t = u = 0, so 



A(s) = 5 l3 8 



u (F{p) , F 2 {p) 



U 4m 2 - m 2 



(C17) 



Thus A(t, s, u) — A(u, t, s) — A(0) and our crossing relations (Eq.(Cll)) give, 

A° = 3A(s) + 2A(0) 
Al = 2A(s) 



(C18) 



from which we have, 



^ = ^(s) + 2A(0)) 



(C19) 



Work presented formerly as part of A.S.J.'s Ph.D. thesis submitted to the Faculty of the Colorado School of Mines in 
April, 1997. 
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TABLE I. Results for numerical solution of cr-model RG flow equations with 6 quark flavors for three different values of the 
UV quark-meson coupling go- The and || indicate that m° q and (fio, Ao) were used to fix m q and (m^,/^) respectively. 



Parameter 


Fit 1 


Fit 2 


Fit 3 


a 


2.967 


3.358 


3.893 


9o 


2.500 


2.765 


3.100 


/xo(MeV) 11 


666.0 


739.4 


818.0 


rn^MeV) 1 


6.54 


6.42 


6.36 


Ao(MeV) 11 


950.0 


937.9 


927.0 


(MeV) 11 


140 


140 


140 


U (MeV) 11 


92.60 


92.49 


92.44 


m^MeV) 1 


7.48 


7.51 


7.54 


m CT (MeV) 
Ai_j,t 
A 3 - pi (100 MeV) 
< > 1 / 3 (MeV) 
alim- 1 ) 

(2a»-5ag)(m- 1 ) 


507.3 
17.3 
43.0 
190 
0.232 
-0.042 
0.677 


536.1 
20.8 
50.8 
195 
0.225 
-0.043 
0.664 


550.4 
30.1 
64.2 
201 
0.220 
-0.043 
0.656 



TABLE II. Comparison of s-wave nn scattering lengths obtained by measurement and various calculations in dimensions of 
inverse pion mass. Experimental results in the first row are from Ref. (34j. The calculations in the second through the third 
row are quoted from Refs. [^, ^5|j3rj| respectively and were performed using the perturbative cr-model, chiral perturbation theory 
(xPT) , and lattice QCD respectively. In the last three rows are the results from our model for the three fits used in Table |l[ 
Comprehensive reviews of nn scattering are given in Ref. section VI-4 and Ref. H. 









2al - 5ag 


Experiment 


0.26 ±0.05 


-0.028 ±0.012 


0.66 ±0.12 


Pert, cr-model 


0.16 


-0.045 


0.56 


X PT 


0.20 


-0.042 


0.65 


Lattice QCD 


0.22 


-0.042 


0.65 


Fit 1 


0.232 


-0.42 


0.677 


Fit 2 


0.225 


-0.043 


0.664 


Fit 3 


0.220 


-0.043 


0.656 
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FIG. 1. Tree diagram 




p(100MeV) 

FIG. 2. RG boson potentials for the a-model. The parameters of the calculation are displayed in Table 3.1 
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FIG. 3. Quark mass funcions for the a- model. The parameters of the calculation are displayed in Table 3.1 
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